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1. Introduction 

A massless quantum field theory has no gap in the excitation spectrum. This can 
be seen, for example, in the Laplace representation of Green's functions. This of course 
does not imply scale invariance; in general properties will interpolate between those of the 
two different conformal field theories describing the UV and IR fixed points. A standard 
(although a little marginal) example of such theory is the 0(3) non-linear sigma model 
with topological angle O = tt, which has central charges cuv = 2 and cm = 1. 

In 1+1 dimensions, many massless theories are integrable. Such theories include well- 
known statistical-mechanical models like the continuum limit of the XXZ spin chain and 
the Kondo problem. Many more are provided by appropriate perturbations of conformal 
field theories. Their study is interesting for several reasons. A few properties are acces- 
sible experimentally; see for example Features of academic interest include Green's 
functions with different anomalous dimensions in the UV and IR, the consequences for the 
topology of the space of relativistic quantum field theories, a better understanding of the 
second law of thermodynamics associated with renormalization group trajectories, and a 
way of understanding perturbations of IR fixed points by irrelevant operators. 

In these lectures we will discuss the scattering theories associated with integrable 
massless quantum field theories. In a massless theory the excitations should consist of 
right-moving and left-moving particles with p = ±E, where we set the speed of light to 
be 1. S'-matrices describing the "scattering" of such particles were calculated long ago 
in for the XXX model. Such objects do not make much sense in traditional S'-matrix 
theory where one requires the existence of in and out states; it is difficult for instance to 
imagine a physical process that would lead to scattering between two particles moving in 
the same direction at the speed of light. Massless S'-matrices in 1 -|- 1 dimensions seem 
to make sense only in the context of integrable quantum field theories. In this case the 
scattering is completely elastic: momenta are conserved individually. We build states by 
acting with creation and annihilation operators on the ground state. These operators have 
non-trivial commutation properties (the Zamolodchikov-Faddeev algebra @,^) encoded in 
the ^-matrix (suppressing internal indices describing the particles): 

n+{9i)n+{92) = si9i - 92)n+{92)n+i9i), (i.i) 

where by convention lZ'^{9i)lZ'^{9j) creates a plane wave with Xi < Xj. This formal 
definition of 5", directly inspired by the Bethe ansatz equations, makes sense in both 
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massive and massless cases. Another way of describing this definition is as a matching 
condition on two-particle wavefunctions. 

Some care must be taken when defining massless integrable theories. For instance, the 
usual proof that an infinite number of conserved quantities implies factorized scattering 
relies on the possibility of separating wave packets in general Q, which is not possible in 
the massless case. Analyticity properties are also not completely clear. Since the particles 
are massless, they are either right- or left-moving. Because S'-matrices for pure left-left 
or right-right scattering can be obtained by a limiting process from physical ^'-matrices 
acting on massive particles, one requires from Sll and Srr exactly the same properties 
as for physical S'-matrices. As we will discuss, the case of Slr is subtler. In general 
these properties can be derived using the definition based on Bethe ansatz wavefunctions, 
without reference to any in and out states. 

After having given the caveats, we would like to explain why finding massless S'- 
matrices is a worthy endeavor. Even if the 5'-matrix found has no physical meaning 
in its own right, many quantities calculated from it do. In these notes we will show 
how to calculate the free energy at non-zero temperature using the thermodynamic Bethe 
ansatz (TBA). Modular invariance relates this to the Casimir energy on a circle, giving 
a "c-function" , which for unitary models shows the evolution of the number of degrees of 
freedom in the fiow from ultraviolet to infrared. This calculation can also be modified to 
give some excited-state energies as well, which give conformal dimensions in the critical 
limit. We will also show how to obtain the ground-state energy at zero temperature in a 
background field, a result related to the chiral U{1) anomaly in the critical limit. Finally, 
we present here some new results on obtaining higher-spin Virasoro conserved charges from 
the massless scattering. 

These results in fact suggest that some aspects of a conformal field theory can be 
described by a theory of massless particles with no left-right scattering This can be 
seen as follows. In these massless but not scale-invariant theories, we have a mass scale M; 
M = gives the ultraviolet fixed point while M oo gives the infared one. The momenta 
and energy of the particles are parametrized by 



E = -p = 



E =p 




for right movers 



for left movers 
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A Lorentz-invariant 5"- matrix element Sll describing scattering of two left movers depends 
only on the ratio of the two momenta, so it depends only on 61—62 and not on M. The left- 
right scattering also depends only on the rapidity difference but does depend on M, because 
the only Lorentz invariant is s = (pi +P2)^- We can always rescale M ^ 00 by shifting the 
rapidities. The LL and RR 5'-matrices are independent of this shift (although Slr is not), 
so they are characterized solely by properties of the infrared fixed point. In this sense one 
can think of the LL and RR S'-matrices as being the S'-matrices for the conformal field 
theory. This should not seem bizarre — many properties of four-dimensional field theories 
(even ones with massless particles like QED) are described by particle theories! 

Another reason for studying massless S'-matrices is that finding them is often an easier 
task than doing the full Bethe ansatz, a result of the constraints of an integrable theory. 
In the massive case, this has become a highly-developed art (see P|-|TI1[| for reviews), and 
many of these lessons can be applied to the massless case. The basic method is to guess 
the particle content based on the knowledge of the symmetries of the problem (and on 
the Lagrangian, if one is known), and then find the simplest S- matrix consistent with 
these symmetries as well as the criteria of factorizability, unitarity and crossing symmetry. 
In many cases, such an S-matrix is the correct one, as can be checked by a variety of 
methods. The symmetries, in particular affine quantum group symmetries, must however 
be analyzed carefully. 

In addition to the XXX model, massless S-matrices have been found for a number of 
models. Continuum theories include the flow from the tricritical Ising model to the Ising 



model [|TT1, the 0(3) sigma model at O = tt and the SU{2)i principal chiral model 0, 
the flows between the minimal models |T^, the Kondo problem the "sausage" sigma 



model [14], and the Landau- Ginzburg flows to the N=2 minimal models 15 . Lattice 



models include integrable higher-spin XXX chains fT^ and the Hubbard model |[T3,|T8| . 

The purpose of these lectures is to provide a pedagogical introduction, so we will skip 
many technical details. We have tried to make the sections reasonably independent of 
one another so that they can be read separately. In Section 2 we give a simple example 
of a massless fleld theory, the sine-Gordon model with imaginary potential. Section 3 
contains a discussion of the usual Thirring model (for properties we discuss, Thirring and 
sine-Gordon can be used interchangeably) and its massless limit. Its main purpose is to 
introduce the physically odd idea of left-left or right-right scattering between massless 
particles from the Bethe ansatz view point. Section 4 gives a simple introduction to 
the thermodynamic Bethe ansatz carried out with the example of massless scattering. 
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Section 5 explores a little more how one can describe certain aspects of conformal field 
theories by massless particles with no left-right scattering. Previously unpublished results 
regarding the non-perturbative analysis of conserved quantities in conformal field theories 
are presented. Section 6 contains some comments on integrable theories with nontrivial 
left-right scattering. Section 7 is based on [0. There we carry out explicitly the study 
of the sine-Gordon model with imaginary coupling and background field, the latter being 
introduced to get a simpler calculation. We show that even if massless scattering appears 
a little odd physically, it at least provides the proper analytic continuation of physical 
quantities in appropriate directions of the parameter space. Section 8 contains conclusions 
and questions of interest. 



2. A simple example of massless field theory 



A simple but generic example of integrable massless field theory is provided by the 
sine-Gordon model with imaginary "coupling" (i.e. prefactor of the cosine term). Recall 
the hamiltonian 



H 



ds 



^U^ + ^{dcl))^ + X cos (3sg4> 



(2.1) 



where we set fHfjy = 8tvj^ at the UV fixed point. This model is well-studied for A real 
and known to be a massive field theory with a trivial fixed point in the IR. Taking A 
imaginary does not look too physical at first. There are however several reasons for doing 
so. The massless fiows between minimal models are obtained as reductions of this model, 



and are unitary even though ( |2.1| ) is not. In addition, (|2.1|) describes the fiow of the 



0{n) model to its low-temperature phase [jT^; this covers interesting physical situations 
in condensed-matter physics such as self-avoiding polymers. 

To see the difference in behavior between A real and imaginary consider the large-t 
limit, where the cosine term is almost marginal and reliable perturbation theory at order 
1/t can be carried out following the analysis of the XY model PDI. The RG equations are 



and 



dX _ 
dp ~ 

1 d(3^ 



47r 



A, 



-TT A . 



(2.2) 



(2.3) 



Prom these one deduces, at first nontrivial order 



db 



4tv 



8n 



1 

t2 



(2.4) 



For A real, the initial derivative of (3 is negative, and the coupling fiows to zero at large 
distance. For A imaginary, the initial derivative is positive and the coupling increases 
monotonically to the IR fixed point with 



0: 



IR 



8n- 



t-1 



t 



(2.5) 



Correspondingly both the UV and IR fixed points are Gaussian models with different radii 
of compactification, and we have a fiow "within" c = 1. 

The most interesting aspect concerns the evolution of the running central charge. The 
latter can be defined in several ways away from the fixed points, for instance using the two- 



point function of the stress energy tensor [^T| or finite-size effects [p2|p3| . Qualitatively, 
these functions should behave in a similar way and are usually expected to describe the 
evolution of the number of degrees of freedom. By analogy with the second law of ther- 
modynamics we expect that such a function should decrease when following a RG fiow. It 
has been proven that the first type of c-function always decreases in unitary theories (the 
c-theorem |^1|), and all known unitary examples of the second type also decrease. In our 
non-unitary problem, it is easy to compute them at first non trivial order in 1/t where 



they coincide. One finds then [|12 



c= 1 + 



12 e4(&-6o)A(^i _ e4(&-''o)/t) 



(2.6) 



^3 (1 +e4(6-bo)/t)3 

that has a roaming behavior. It indeed has c = 1 in the UV and IR but goes up and down 
in between, reaching a pair of extrema with values 

2 1 



c± = 1 ± —p^ ■ 



(2.7) 



Of course the fact that c can increase and does not obey a nonunitary version of the c- 
theorem is not surprising and is an obvious consequence of the imaginary coupling. One 
might hope to slightly modify the perturbation to obtain an exact fiow which would stop 
at the value c_|_. This does not seem impossible in view of [^, which describes fiows that 
exactly interpolate between minimal models and flows that go "very near" them. 

The fact that c is not monotonic sheds some doubt on the reliability of the running 
c-function as a measure of the number of degrees of freedom. Let us emphasize that ( |2.6| ) 
is indeed related to the absolute ground state of the theory — it is the ''Cef/' supposed to 
qualitatively replace c in the non-unitary case. 
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3. Massless scattering at the conformal point 

We consider the example of the massive Thirring model on a circle of length L with 
hamiltonian 

H = j dx [-1 {tpid^tpi - il^t^xil^'i) + rriQ {tp^ij2 + ^^V'l) + '^9o'^ttpt'^2''Pi] ■ (3.1) 



which is solvable by the Bethe ansatz [p5l , |26[| . The analogous calculation is done for the 



sine-Gordon model in |^^. As a result eigenenergies take the form 

£^ = mo cosh^i, (3.2) 

i 

and momenta 

P = mp sinh^t, (3-3) 

i 

where the are bare rapidities of pseudoparticles satisfying the Bethe ansatz equations 

exp(zmoLsmh6) H sinh(6 - - 2z^) = ^' ^'"'^ 

and cot/U = —^Qo- One recognizes in (|3.4| ) the conditions for the wave function to be 
periodic, the phase shift being a combination of a free term and factorized scattering 
between pairs of pseudoparticles. Thus the elements of the product in ( |3.4| ) are the "bare" 
iS-matrix elements for the pseudoparticles, which we denote as So{^i — ^j) = exp[z(/)o(^i — 
^j)] = jz^- Equation (|3.4| ) follows from the form of Bethe wave functions 

V^(xi,...,a;Ar|^i...,^Af) = exp(imo^XiSinh^i) + zA(^i -^j)sign(^i - Q]- (3-5) 

i i<j 

This wave function is almost free, with only a phase shift when the coordinates are ex- 
changed. Of course since we are dealing with fermions, the real wave function has still to 
be antisymmetrized so 

^'(xi, . . . , xn) = ^ sign{P)tp{xi, . . . , xn\^pi, ^pn), (3.6) 
p 

(this makes sense because A is an odd function, or equivalently Sq{^)So(—^) = 1) and all 
rapidities must therefore be different. 
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The study of the quantum field theory requires building the ground state by filling 
the appropriate Dirac sea, and then finding the excitations and their scattering. We must 
mention a difficulty concerning the choice of UV cutoff. In the attractive regime go > 
essentially all cutoffs produce the same results, but the situation is different in the repulsive 
regime go < 0. (This is unfortunately the regime where one truncates the sine-Gordon 
model to describe massive perturbations of Virasoro minimal models by the (pis operator 
||28|| .) The rapidity cutoff of |2^] leads to many more particles in the spectrum than the 



cutoff of [pO[, which is essentially a lattice cutoff [31]. The latter regularization reproduces 



the conjectures of M and seems in agreement with the expected physics of perturbed 



minimal models [^2|. We avoid entering into technical details and give a "morally" correct 
discussion in the following. 

Another difficulty arises in the eigenfunctions of as a result of terms like (sign S). 
Depending on the regularization, different relations between the bare coupling go and the 
parameter /i are found. In ||2^] the relation /i = is used instead. We shall carry out 



the discussion in terms of the variable fi. 

The ground state is easily built by filling the sea with ant ipseudopar tides, i.e. filling 
the line /m(^) = tc. There are various kind of excitations. We shall discuss only the 
case of solitons s and antisolitons a (these are the only particles with non-vanishing U{1) 
charge). For instance a pair of solitons is simply obtained by making two holes in the 
ground state distribution. Antisolitons, or pairs of solitons and antisolitons are obtained 
similarly, with the addition of some strings of pseudoparticles around /m(^) = 0. Of 



course the introduction of such holes induces a shift in the distribution of the ^'s, the 
so-called backflow. As a result the mass and rapidity of the solitons are renormalized, 
giving 

i? = mcosh^, p = msinh^, (3-7) 



with 



^e^(i-7) 

m = mo—, -tan7r7, 6* = 7A, (3. 

7r(7- 1) 



with 7 = ^ and A denotes the bare rapidity of the particles (s, a) (again there are some 
slight differences between authors for this formula). As shown in the S'-matrix of 
solitons can be extracted from the Bethe ansatz equations. The 5'-matrix elements are 
defined as follows. Suppose first we have only one soliton, with bare rapidity Ai. Then 



the total phase shift coUected by the wave function when the argument of the sohton goes 
around the circle is 

(j)! = mpLsinhAi + (pojXi - Cj) (3-9) 

j 

where the sum is taken over all pseudoparticles in the sea and ^ indicates shifted (with 
respect to the ground state) rapidities due to the backflow. Suppose then we have two 
solitons with bare rapidities Ai and A2. Then the total phase shift of the wave function 
when the argument of the first soliton again goes around the circle is 



LsinhAi + ^(/)o(Ai -f,). (3.10) 



mo 



where ^ are shifted rapidities. One then defines the S'-matrix element by InS" = i{(p2 — (pi)- 
Complete computation shows that it depends only on the difference of the rapidities. 
Moreover one also checks that for more particles, the phase shifts simply add and the 
scattering can be decomposed as a succession of two-particle ones. The resulting S'-matrix, 
therefore a solution of the Yang-Baxter equation, has the well-known matrix elements 

a = Z{9) sinh (^^^^) , b = Z{9) sinh (^^^ , c = Z{e) sinh (^y^ (3.11) 

where a corresponds to ss ss scattering, b to sa —>■ sa and c to sa — as. We give the 
the expression for normalization factor Z in sect. 7. The symmetry under s ^ a gives the 
remainder of the elements. We have parametrized 

T-f-rfr. (3.12) 

2(U t + 1 



The S'-matrix (|3.11|) can be manipulated to exhibit Uqsl{2) symmetry [33] with 



g = -exp(^-yj. (3.13) 

This is a dynamical symmetry; there is also a kinematical symmetry UqgSl{2) with go = 
— exp^— following from the Bethe ansatz equations [34]. Notice the shift of the 
denominator. 

In order to reach the deep UV limit, we let the mass mo go to zero. Particles with 
non-vanishing energy must have rapidities with very large modulus, of the order ^0 = 
ln(M/mo) >> 1, where M is a not-yet-defined parameter with the dimension of mass. 
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There are thus two regions of interest in which we set respectively ^ = ^o+^ and ^ = —^o+O, 
9 remaining finite. The spectrum obviously splits into right and left excitations with 

En - ^e^"e^ ^ = ^^'^ " " f^"' = "^^^ ^'"''^ 

and we have a doubling of species ai,^/j, s^.i? (see 1). For pseudoparticles of the same kind, 
the phase shifts are unchanged as — = 6i — 9j . For particles of different kinds however 
— ^ ±2^0 ioo so the phase shifts become constants, independent of the rapidities. 
As a result, in the computation of the S'-matrix, the LL and RR scattering are the same 
as the ones for corresponding massive particles computed above Sll = Srr = 5*0, while 
the LR scattering becomes trivial. 

The iS-matrices in the massless limit can also be obtained by studying the XXZ spin 
chain, which is a lattice regularization of the massless Thirring model. This is a simple 
generalization of the work of on the XXX chain. 

It should be clear finally that the properties of a massless theory with trivial LR 
scattering are independent of the mass scale M. Indeed, changing M is equivalent to 
shifting and the analysis only depends on rapidity differences. 

4. Thermodynamic Bethe ansatz 

The technique we now discuss involves computing the free energy of an integrable 
lattice model (e.g. the XXZ model) or an integrable quantum field theory (e.g. the Thirring 
model) on an infinite line at finite temperature T. There are two approaches to this 
calculation. In the traditional "bare" approach, one finds the energy and entropy of the 
states of the model using the Bethe ansatz. The thermodynamic state is the state which 
minimizes the free energy. The limit T ^ gives the ground-state energy and the vicinity 
T ^ the structure of low-lying excitations. In the massive Thirring model this approach 
starts with the bare equations ( |3.4| ). In the second approach one forgets the bare theory 
completely and studies instead the thermodynamics of a gas composed of the various 
"physical" excitations (like the soliton and antisoliton of the sine-Gordon model) scattering 
with their respective 5"- matrices. As in the first approach, one determines their energy and 
entropy and again minimizes the free energy. The first approach was used in works like 

^ In some cases like in there is an additional phase in the definition of Sll and Srr. 
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|3^,0. The second approach, pioneered in is more recent, and is usually called the 
thermodynamic Bethe ansatz (TBA). 

The second approach allows some convenient short cuts: instead of solving the theory 
one first establishes that it is integrable, conjectures the excitations and their S'-matrix 
using intuition and symmetry arguments (with some care), and uses the TBA to derive 
various properties. When both approaches can be implemented, they of course give the 
same results; the quasiparticle excitations used in the second can be found from the Bethe 
ansatz equations by filling the Fermi or Dirac sea. However, there are examples (at least 
in the massive case) where a lattice model is not integrable but its continuum limit is an 
integrable quantum field theory. Usually the only way of defining the continuum model 
is by a perturbed conformal field theory, so the usual Bethe ansatz methods cannot be 
applied; the only recourse is to use the second approach. A classic example is the Ising 
model at T = Tc in a magnetic field ^H. 

As a simple example we describe the TBA for a single type of massless particle, say 
right-moving, with energy and momentum parametrized as in ( |1.2| ). The scattering is 
described by a single S'-matrix element Srr. Quantizing a gas of such particles a circle of 
length L requires the momentum of the zth particle to obey 



exp 



[^^Ljl[SnR{e,-e,) = l. (4.1) 



One can think of this intuitively as bringing the particle around the world through the other 
particles; one obtains a product of two-particle S'-matrix elements because the scattering 
is factorizable. This is the renormalized equivalent of the bare relation ( |3.4| ). 

Going to the L ^ oo limit, we introduce the density of rapidities indeed occupied 
by particles p{9) and the density of holes p. A hole is a state which is allowed by the 
quantization condition (f4.1| ) but which is not occupied, so that the density of possible 
rapidities is p{9) + p{0). Taking the derivative of the log of ( [4.1|) yields 

Mr 

2n[p{e) + m] = -^e' + J ^e- e')p{e'), (4.2) 

where 



^ Observe however that there is another integrable lattice model based on Eg that is integrable 
and has the same scaling limit as Ising in a magnetic field. 
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To determine which fraction of the levels is occupied we do the thermodynamics . The 
energy is 



and the entropy is 



M 

p{e)-e'de, 

/oo 
[{p + p) ln(p + p) - p ln(p) - p ln(p)] dO. 
-oo 



The free energy per unit length T = {£ — TS) /L is found by minimizing it with respect 
to p. The variations of £ and S are 



5£ 
5S 



°° M 
Sp—e'de 
^ 2 



- oo 
oo 



[{6p + dp) ln(/9 + p)- 5pln(/9) - 5p\n{p)] dO. 



It is convenient to parametrize 



P{0) _ ( 6 



(4.3) 



givmg 



5S 



5p In f 1 + e'/^) + 5p In f 1 + e''/ 



d9. 



Using (|4.2|) allows us to find p in terms of p. Denoting convolution by this gives 
27r{5p + Sp) = $ ★ 5p so 



ss 



T 2n 



6pde. 



Hence the extremum of JF occurs for 



—e^ = t + T— ★ In ( 1 + e"'/^ 
2 27r 



(4.4) 



and one has then, expressing p from ( [4.2|) and using ( ^I^ ) 



(4.5) 



The ground state energy £q cannot be obtained by this method since all the information 
we use is the structure of excitations above the ground state, so we set £o = Q for the rest 
of this section. 
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The limit T ^ of this system is interesting. We introduce the positive and negative 
parts of the pseudoenergy satisfying therefore 



M 



2^ 



★ e 



(4.6) 



In this hmit the solution is e~ =0, e"*" = ■ K follows from this and (|4.3|) that p ^ 
as T ^ 0, which is required because our TBA provides the structure of excitations over 
the ground state. In general — e~ (resp. e+) gives the excitation energy for holes (resp. 
particles) . 

The knowledge of JF leads to the determination of the central charge of the theory. 



We have JF 



In Z, where Z is the partition function of the one dimensional quantum 



field theory at temperature T (in the following we refer to this point of view as "thermal"). 
In Euclidean formalism, this corresponds to a theory on a torus with finite size in time 
direction R = 1/T. By modular invariance, identical results should be obtained if one 
quantizes the theory with R as the space coordinate . For large Z = ^-^i^)^ ^ where 
E{R) is the ground-state (Casimir) energy with space a circle of length R. Thus T = 
E{R)/R. In the following we refer to this as the "finite-size" point of view. Conformal 
invariance requires that at a fixed point this Casimir energy is E{R) = — f^, where c is the 
central charge [^. Going back to the thermal point of view, JF = — ""'^^ and the specific 
heat is C = 

Observe from ( [4.4| ) and (j^^ that the free energy does not depend on the mass scale 
M, because it can be rescaled by a shift in rapidities. By dimensional analysis one has 
therefore JF = x constant. This scale invariance is a manifestation of the fact that 
Sll and are describing only conformal properties. With massive particles or with 
nontrivial left-right massless scattering, JF does depend on M/T, giving a running central 
charge. 

We can analytically find this central charge from ( [4.4|) . We take the derivative of 
with respect to 9 and solve for e^. Substituting this in ( [4.5| ), we have 



)/T 



T 


^ de 


'de 










_d9 


T 


r de 




/ dO- 




2^ „ 


1 dO L 




T 


^ de 




2^ J 





jde'Hi + e-^^y-me-o')^^^ 



Infl + e-'/^) + (e- —e^) - 



de 
dO 



ln(H-e-^/^) + 



1 + e^/T 
12 



where we use (^4.4| ) again to get to the second hne. We can replace the integral over 9 with 
one over e, giving an ordinary integral 



^=-7-1 de 

47r -/.(_oc) 



1 + e^(^)/^ 



A change of variables gives 

where L{x) is the Rogers dilogarithm function 



and xq = exp(e(— oo)/T) is obtained from (^]J) as 



Xo \ Xo J 

with / = 7^ f $. 

For example, when the S matrix is a constant, $ = 0, xq = 1 and 

where we used L(l/2) = j^. Here we find cl = \- In a left-right- symmetric quantum field 
theory, the right sector makes the same contribution, giving the total central charge c = | 
required for free fermions. 

For the nonunitary Lee- Yang iS-matrix [38] one has / = — 1. In that case xq = "^^^ ■ 
Using L {^—^^ ~ T,\ ^'^^ finds cr = and c = | after left and right contributions 
have been collected. This is indeed the effective central charge for the Lee- Yang problem; 
it is not equal to the true central charge c = — ^ because of the presence of an operator 
with negative dimension in the vacuum. 

It is possible to do the same computation with the massless pair (a, s) scattering with 
The computation is technically more complicated because the scattering is non- 
diagonal. We just refer the reader to references |]3^,0 for details. Simply observe that in 
the case A* = f the scattering becomes diagonal and because of the doubling of the number 
of species, the above calculation gives rise to c = 1 as expected. 

Besides the central charge, some conformal dimensions can also be identified using the 
TBA. To do so one includes an imaginary chemical potential /i^ (not to be confused with 
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11 in section 3) for each species of particle h. As before, we minimize the corresponding free 
energy Q = £ — TS — X]& A^fo-M?, and find similar results with, in most cases, exp(— e^/T) 
replaced by exp( — (ef, — iit,)/T). In the finite-size point of view, the introduction of a 
chemical potential amounts to considering the theory on a circle of length R with twisted 
boundary conditions. As is well known, the ground state energy in that case gives an 
effective central charge Ce// = c — 24h where h is related to the twist. For RR scattering 

2 

given by ( |3.11j ) for instance, with Hs = —^a = ianT/t one finds h = j^^j^. Thus one 
recovers the dimensions of vertex operators in a Gaussian model. 

The question of reconstructing the whole quantum field theory from a massless scat- 
tering theory with no left-right scattering is still open. For Virasoro minimal models two 
independent such theories are probably necessary, as there are two fundamental quantum 
groups, or two labels in the Kac table. A bit of progress in this direction is presented in 
the next section. 

5. Integrable CFT and massless scattering: Virasoro conserved quantities 

We showed in sect. 4 how one obtains the free energy in a massless integrable theory, 
and compared this result with conformal field theory predictions. In an integrable theory, 
the energy S = (E) is just the first of an infinite series of conserved charges. These 
conserved charges can often be expressed as suitably regularized powers of the energy- 
momentum tensor. Since we derive the particle densities of the thermodynamic state, we 
can calculate the expectation value of any quantity which can be expressed in terms of 
the particles. This suggests that the conserved charges should be related to expectation 
values {E"^). In this section we show that in the sine-Gordon model, this is indeed true, 
thus verifying non-perturbatively what had been shown classically and perturbatively in 
the quantum theory P0| , |il[1 . 

We start with free fermions with antiperiodic boundary conditions on a circle of length 
R in order to select the ground state. We are thinking about the system in the finite-size 
point of view discussed in the last section. Consider the quantity 




(5.1) 
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The sum can be evaluated by C-function regularization leading to 

2fc+l 



(5.2) 



and 



{E^^)gs = 0, (5.3) 



where the last result follows from C(— 2/c) = 0. For k = one gets S = (E) = with 
c = 1/2 (c/2 appears because we concentrate on one chirality). 

We can compute the analogous quantity from the thermal point of view by putting 
the particles on a circle of large length L at temperature T = 1/R. The TB A analysis of 
sect. 4 gives 



{E^TBA = j de (^^j p{6) 



oo / ^ ^ f)\ n 



oo 

■oo 



where we used the fact that 27r(p + p) = LTde/dO, which is proven by showing that they 
obey the same integral equation. 

For free fermions, e = Me^/2. The expectation value ( p.4|) is 

T 

^ ^ n / n-\ 



{E^)tba = -BirR'"" / ln(l + e-^)dx, 

Jri Ztt 







or, restricting to n = 2k + 1 



{E''^')tba = ^^-^h±niR-2i^-^ _ _L_^ c(2fc + 2). (5.5) 



R 2 



To compare ( |57^ ) and ( |5.5| ) recall the identities ^ 



where are Bernouilli numbers. Hence we find 

{E''+')TBA = i-l)''+'^{E''+'),s. (5.6) 

Jri 
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For /c = we recover T as derived in the previous section, using the fact that E = —T 
for theories with no left-right scattering. More generaUy (|5.6|) follows from the relation 
of {E'^^^^)tba (resp. {E"^^^^) gg) to the energy-momentum tensor component T^^^ (resp. 
Tyy^) and T^x + Tyy = at the conformal point. 

Notice that when n is even, the results are quite different, since {E^'^)tba 7^ while 
( p75| ) holds. The usual argument is that such even powers do not correspond to any local 
quantity in the quantum field theory, and therefore the two results cannot be compared as 
we did in ( p.6|) . 

The generalization of this computation to the case of nontrivial scattering is not 
straightforward, but it is useful. We can compute - at least numerically - the quantities 
{E'^^^^)tba from (|5.4| ). In general the conformal field theory cannot easily be described by 
oscillators as in the free theory above, so we do not compute the equivalent of {E^'^'^'^) gs. 
What we can however compute using only the Virasoro algebra are quantities like 



where the integral is over the period of the cylinder, and average is taken in a state of 
conformal weights [h, h), generalizing the ground state. The double dots indicate normal 
ordering on the cylinder. Recall that this normal ordering leaves room for non-vanishing 
expectation values. These can be computed by explicitly performing the subtraction of 
the divergent terms. The result coincides with the simpler zeta regularization. Such a 
quantity cannot be directly compared to {E'^^^^)tba because of a non-trivial renormal- 
ization factor. In the free fermion case for instance (J : :)gs = {^)^ ||C(~3) while 



{E^)gs = {^f ^C(-3), because : := | : d^tPdtP : - ^ : d^tPdtP :. However this nor- 



malization factor occurs from short-distance singularities and therefore does not depend 
on boundary conditions. We can therefore compare ratios of the moments of T, dT, . . . 
and {E"')^s for different boundary conditions around the circle, or equivalently the choice 
of state in which the average ( |5.7| ) is taken. As explained in the previous section, in the 
thermodynamics this means taking imaginary chemical potentials for the particles. 

A crucial point is that we must treat the theory as a non-minimal, non-unitary theory 
with central charge c = 1 — jj-^ii^- We use this value of c, denoted as Cmin, in the Virasoro 
algebra computations. Boundary conditions other than these twisted ones give rise to 
Cef f = Cjnin — 24/i, where the conformal dimension h is computed with respect to the Cmin 
ground state. The usual sine-Gordon ground state (with Ce// = 1) is then interpreted 




(5.7) 
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as arising from an operator of negative dimension. With these definitions one finds for 
instance, using the Virasoro algebra and ^-function regularization, 

{j ■.T':)h = ^{10clff + 40c,ff+4crmn), (5.8) 

, /■ o , 1 s 2 192 32 6 \ , , 

{J - T •.)h = -J2^ y^eff + '^'^Ceff + — Ce// + —Cmin + ^C^mCe// 1 , (5.9) 

( j : Td'T = (^c^// + f c^..) . (5.10) 

The Ce// result from the non-zero {Lo)h on the cylinder. 

Our strategy has been simply to compute numerically the values of {E'^^^^)tba for 
various chemical potentials, denoting these by {E'^^~^^)of For sine-Gordon the chemical 

2 

potentials are Hs = —^a = ianT/t; this results in Ce// = 1 — 24/i, where h = j^^^- 
The general TBA equations with fugacities are written out in [^. Since the numerics 
are crucial to obtaining our result, we describe the methods briefly. The multi-function 
generalization of ( [4.4| ) is of the form 

ea{9) = MO) + d0'(pab{9 - 9') ln(l + Xbe-'"^''^/^) 

b 

To flnd ea numerically, we solve this iteratively. We guess the initial e^; = fa usually 
works. Then one evaluates the right-hand side numerically by discretizing the integal; this 
gives the next guess for e^. Usually the iteration converges; occasionally one needs to use 
take a linear combination of a;(guess)-|-(l — a;) (iteration) for the next guess. More elaborate 
methods to improve convergence are described in Once this procedure has obtained 
ea to the desired accuracy, the expression ( p^) for {E'^) can then be numerically evaluated. 
We note that this numerical procedure for solving non-linear integral equations is generally 
far simpler to implement that those for solving non-linear differential equations. 

The numerical results are rather interesting. As before, we have Cmin = 1 — 6/t(t + 1) 
and Ce// = 1 — 6a^/t(t -|- 1). We flnd 

(E^)^ =ft{104ff + 40Ce// + 4c^^n) 

16 8 (5-11) 

=gt{clff + I2clff + 40Ce// + 2Ce//C^in + yC^*n + ^C^n)' 

Moreover, we flnd that at least for {E^), the prefactor takes reasonably simple values; to 
excellent numerical accuracy we have 

7 1001 ^ 27r^ 

- 48 70 ~ 7200 ~ 143 
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This is a hint that these numbers can be derived analyticaUy from the TBA, but we have 
tried and failed to do so. We can make the amusing observation that for free fermions 
can be written in terms of polylogarithms, just hke the dilogarithms written in the last 
section for {E). We also note that the free energy of the impurity in the Kondo problem 
(see IjlJ) can be written as a sum of these conserved quantities, another hint of interesting 
hidden structure. 

Hence we can fit our numerical results to ( |5.8|) for n = 3 and to a linear combination 
of ( p.9| ) and (|5.10| ) for n = 5. We see that (|5.8| ) is indeed proportional to {E^)tba and that 



{E')^ oc (y : : + j^^^!!!!^ j Td'T :),. (5.12) 

After the numerical computation was completed we checked that ( |5.12| ) agrees with the 
conserved quantity at grade 5 in [Q. This is of course no surprise. Recall that in the 
classical sine-Gordon theory, the conserved quantities are precisely expressed as the sum 
of odd powers of the momenta: we simply check here that this result holds in the quantum 
theory as well. This is expected, but as far as we know, was checked only perturbatively 
so far 1^,^. Hence by massless scattering we recover not only the central charge and 



conformal weights of a conformal field theory, but also the conserved quantities which 
involve the Virasoro algebra itself, making the connection between the two points of view 
a little closer. One might wonder if there is an action of the Virasoro algebra on the 
massless particles. 



6. Nontrivial left-right scattering 

The massless Thirring model has no non-trivial LR scattering because L and R excita- 
tions are infinitely separated in the rapidity plane. On the other hand, the most interesting 
situations occur when the LR scattering is non-trivial. In that case, the theory is not scale 
invariant, and is described by two different conformal field theories in the UV and IR 
limits. To get such a situation in the Thirring model, we need to arrange for massless L 
and R excitations that both occur around the same region of rapidities. A way to do so is 
to choose a purely imaginary hwe mass in ( |3.1| ) mo = —i\mo\. Indeed the result ( p7^ ) still 
holds, so 

E' = — z|mo| ^^cosh^i. (6.1) 

i 
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If we restrict to the consideration of real energies we need Im{^) = ±n/2. Then for 
^ = ±i|^ + z^, e = ±|mo| sinhz^. One therefore expects the ground state to resemble figure 
2 with the half lines 

i?e(0 < 0, ImiO = |; i?e(0 > 0, lm{0 = -|, (6.2) 

filled up. Actually, because the theory is not free, the determination of the ground state 
is slightly more delicate — the interaction between the various Bethe ansatz roots must 
be considered. The choice of the cutoff is also important, as well as the sign of qq. One 
finds typically that the picture (|6.2|) is almost correct, up to some exponentially decaying 
density on the other side of the half-lines. This produces therefore the necessary massless 
excitations around a common region ^ = 0. More details of this approach will be presented 
in 



m 



An imaginary mass in the Thirring model is like imaginary prefactor in front of the 
cosine term in sine-Gordon, so we recover the situation discussed in sect. 2. As explained 
there the appearance of imaginary numbers is more natural than may appear at first 
sight. Just as the massive minimal models perturbed by (pis are related to the ordinary 



sine-Gordon model |2S], the massless fiow between minimal models [21,23] is related to 
the model with imaginary mass. The non-unitarity of the imaginary-mass model does not 
exclude unitarity for a subsector (the perturbed massless minimal model) . For more details 
see [0. 

It is easy to generalize the TBA of section 4 to models with a a non-trivial Slr- This 
time, the running central charge depends nontrivially on M/T. Its UV and IR values can 
be easily found. As discussed in the introduction, we expect that the IR conformal field 
theory is characterized by only Sll and Srh, so its LR scattering should be trivial. One 
finds as before 



Cm = cr + cl = 2cr = -^L I — ) , (6.3) 



where yo is the solution of 1/j/o = (1 + ^/Uo)^^ with Ii = ^ J ^ll- In the UV coupling 
between left and right particles has to be considered leading to 



6 

cuv — 



2L 



1 + xi J V 1 -t- xo 



1 



(6.4) 



where 1/xi = (1 + l/xi)^i+-^2 ^^d I2 = ^ J ^lr- 

A simple example of a massless integrable field theory is the flow from tricritical to 



critical Ising model. As discussed in jlll] the spectrum consists of a right mover and a left 
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mover, the Goldstino resulting from spontaneously-broken supersymmetry. Because the 
IR conformal field theory is a free fermion, Sll and must be trivial. The left-right 
scattering is given by 

SRL{eR - Ol) = -tanh (^^y^ - • (6.5) 

The compatibility of left-right and right-left interchange of arguments of the wavefunction 
requires that 

Slr{0l - Or)Srl{Or - Ol) = 1, (6.6) 

so here 

SLR{eL - Or) = tanh (^^^^ - • (6.7) 

In the IR limit (the Ising model) where 0r — 9l —>■ oo one checks that both matrix elements 
go to —1 as expected. With this ^S-matrix we have /i = and /2 = 1/2, so xq = 1 and 
xi = ^2^^ ■ Using values of dilogarithms given in sect. 4 we find cir = 1/2 and cuv = 7/10 
as desired. 

Following (|6.5|) and (|6.6| ) notice that 

SRL{0)SRL{-e) = -1, (6.8) 

a result that must be carefully compared to the usual S{9)S{—9) = 1 for diagonal massive 
(or left-left or right-right) scattering. 



7. Sine-Gordon model in a background field 

In this section we discuss the sine-Gordon model in a background field coupled to 
the U{1) soliton-number charge. In the traditional bare approach, this field would modify 
the Dirac or Fermi sea. In our approach, this makes it energetically favorable for physical 
particles to appear in the vacuum, even at zero temperature. In the sine-Gordon case with 
positive background field, only the negatively-charged particles appear in the vacuum. 
Their mutual scattering is diagonal, so the problem is technically easier than the finite- 
temperature TBA problem, where both kinds of particles appear in the thermodynamic 
state. We discuss both the cases A real and A imaginary, hence giving a (partially) non- 
perturbative treatment of the problem raised in sect. 2. For a more complete study see 



20 



The sine-Gordon Hamiltonian with a constant external U{1) gauge field is 



H = I dx 

where 



QA, (7.1) 



In the ordinary case A real, Q is the integer- valued soliton topological charge, normalized 
so that the soliton (antisoliton) has Q = 1 (—1). A is a constant with the dimension of a 
length. We then consider the corresponding specific vacuum energy S{A, A) as a function 
of A. As before, we parametrize (3gQ = Sn-j^. 

Before we turn to the scattering theory, it is worth looking at the action ( [7.1| ) from the 
perturbative (in A) point of view. Dimensional arguments as well as explicit perturbative 
calculations show that the background field works as an infrared cutoff at scales ~ A and 
therefore if A ^ ;\(i+*)/2 ^j^g theory is in the ultraviolet regime. As a leading A ^ oo 
approximation we set A = in ( [7.1| ). This theory is the continuum limit of the XXZ model 
in a magnetic field, which has been studied in refs. [^7| - ^ . As can easily be inferred from 
the action, it is a Gaussian model whose radius of compactification depends on A (amusing 
finite-size corrections occur in the related XXZ model in a field due to commensurability 
problems between R and the scale A). Redefining d(f by a shift gives for the ground-state 
energy density 

£(A,0) = -^A\ (7.3) 

At any critical point, S{A,0) is proportional to A^, since there is no other scale in the 
problem. The coefficient is proportional to the chiral anomaly (found from the JlJl 
OPE) [0. 

For A 7^ the scaling argument shows that S{A, A) is a function of the dimensionless 
variable 

e = AM2/(i+*) 
and by parity has a perturbative expansion in 

^pert(6= (7.4) 

TT < * 



TT 

1=0 
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One can use perturbed conformal field theory to derive 

t 



ko 



1 + t 



k2 



4 \ 1+t J ^ ( 2t 



1 + t 



We expect the series ( [7 .41 ) to have some finite radius of convergence ^q, defining therefore 
an analytic function SpertiO 1^1 < ^o- The perturbation theory is the same for A real 
or imaginary, so (|7.4| ) holds all around ^ = 0. However, the scattering theory depends 
crucially on the nature of A: for real A the particles are massive, while for A imaginary 
they are massless. 

Consider first the unitary massive sine-Gordon model (A real in ( |7.1D ) and let m be the 
mass of the corresponding charged particle (soliton). As usual the on-mass-shell momenta 
[E, p) are parameterized in terms of rapidity 9 

E = m cosh 9 ; p = m sinh 9 

In the field ( [7 .11) every soliton (antisoliton) acquires additional energy —A (A). It is 
clear that if A ^ m the state without particles is no longer the ground state. The true 
vacuum contains a sea of positively-charged solitons which fill all possible states inside 
some (A-dependent) "Fermi interval" ~B < 9 < B. The non-trivial scattering of the 
solitons certainly infiuences the structure of the ground-state sea. However, only one kind 
of particle is in the sea (this can be checked more completely jH^,^); for A > this is 



the soliton. The solitons scatter diagonally among themselves, the two-particle amplitude 
being a{9) from ( |3.11| ). 



As in the finite-temperature case, we define the density of particles p and density of 
states p + p. To obtain the ground state energy we minimize 

giRe) _ g{Re) = cosh 9 - A)p{9)d9, (7.5) 

(the superscript Re is added to stress that currently we address the ordinary sine-Gordon 
model with real coupling A) subject to the quantization 

2tt[p + p] = mcosh9 + ^ p (7.6) 
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where the kernel $ foUows from the sohton-sohton S'-matrix and reads exphcitly 

HO) I /■ e^-^sinh^(^^ da; 

loga(^) = / ^:Jntu. 7TZ- (7-7) 



2tt 2mde ^ ^ ' J 2 cosh ^ sinh ^ 27r ' 



The equations (|7.5|) and ( [7 .61) and the equations for B (given by minimizing the energy 
with respect to B) can be put in a more convenient form by defining the function f{9) as 

f{e) = A-mcoshe+ [ de'^ie -e')f{e'), (7.8) 

J-B 



where this equation is good only for \9\ < B. Replacing A — mcosh.9 in (|7.5| ) with this 
and using ( |7.6D one finds that 

/B 
dO cosh^ f{e). (7.9) 
-B 

The boundary conditions f{±B) = determine B. 

We can understand the meaning of the function / as follows. Define e"*" as the energy 
of particle excitations above the ground state, and e~ as the energy of holes. By this 
definition, e"*" > and e~ < 0. A variation of the energy is thus 

S£(R'^)(A) = J (mcosh^ - A)5p{6)d6 = J e+5p - e'Sp. (7.10) 



Using (|7.6|) to reexpress 5p as a function of 5p we find 



5S^R-){^A)= I (e+ + e--^*e-)5p (7.11) 



so by comparing ( |7.10| ) and (|7.11[ ) we find 



$ 

mcosh6' - A = e+ + e" -k e~ . (7.12) 

2ti 

Using this in ( [f?5|) gives 

SiR-){^A) - ^(■^^)(0) = — [ de [cosh^e- + e+p - e' p] . (7.13) 

2% J 

In the ground state p{9) = when e~^{9) > and p{9) = when e~{9) < 0. This means 
that the last two terms vanish, and we obtain (|7.9| ), where / = —e~. 

From a formal point of view one may interpret the whole scattering approach and the 
resulting system (|7.7| )-(fr9|) as a way of summing up the perturbative expansion (|7.4| ) at 
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real ^ and going beyond the radius of convergence along A real. One would like a similar 
tool to sum up (|7]^) at A purely imaginary. To do so we make a guess inspired by the study 
of the massive Thirring model and section 6. We assume that in the sine-Gordon model 
with imaginary coupling or Thirring model with imaginary mass, the number of species 
doubles, so now we have a pair of left and right massless particles. We also assume that 
the LL (identical to RR) and the LR scattering are nontrivial. The energy spectrum ( |1.2| ) 
is gapless. Turn on the background field A > as before. The positively charged particles 
are always excited in the ground state; we assume as in the massive case that they are 
the only particles contributing to the thermodynamics. Now the right- and left-movers 
fill respectively the semi-infinite Fermi intervals —oo<9<B and —B < 9 < oo with 
some Fermi boundary B ~ logA/M. Again it is a straightforward Bethe ansatz exercise 
to derive the following system of integral equations 



B oo 

QA~^ = fR{9)- [ ^LL{9-9')fn{9')^- [ <^rl{9 - 9')pL{9')f^; 

2 J ZTT J ZTT 



-oo -B 
oo B 



(7.14) 



QA-'^ = fL{9)- / $ll(^-^')/l(^')^- / ^RL{9-9')fR{9')'^. 



The positive functions fR{9) and /l(^) are defined in the Fermi intervals oo < 9 < B and 
—B<9<oo respectively, and are restricted by the boundary conditions 

fR{B) = fL{-B) = (7.15) 

Finally, the ground state energy, which we now call £^^'^\A), is evaluated as follows 

(A) (0) = -— / fR{9)d9, (7.16) 

where we have taken into account the obvious symmetry fR{9) = /l(— ^). 

In eqn. ( [7.14| ) we introduced the U{1) charges ±Q of the massless particles, which, 
with the knowledge assumed in these lectures, we cannot fix in advance. It is possible to 
carry out the computation with undetermined kernels, and fix them at the end by requiring 
the result to provide analytic continuation of the perturbative series to A imaginary. To 
save time, let us give the answer and justify it. One has for $ll(^) the same expression as 
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in the massive case, but with a shift t ^ t — 1. For ^rl we have the same, with a further 
shift e ^ e + i^{t -1). This results in 



1 , f e-^inh^^(^ duj 



27T 2m d9 " ^ y 2cosh^sinh^i^27r 
27r 2ni d9 °S ""^^^ ^ j 2 cosh ^ sinh ^^^^ 27r ' 



(7.17) 



Although this new BA system ( |7.14D -( |7.16| ) has a rather different form from that of 
eqns. ( |7.7[ )-(fO|), it is easy to relate the two in the UV region A ^ 00 where in both 
systems B ^ 00. In this limit the right and left Fermi intervals have a broad overlap at 
—B < 6 < B. Near say the right Fermi boundary 6 ^ B (where the main contribution 
to ( [7.16| ) comes from) we can forget about the left one and solve ( |7.14| ) for /l(6') by the 



Fourier transform with B — > 00. The resulting equation for fniO) is 

rQA-^ = fn{9)-j ^9 - 9')fn{e')^ (7.18) 



where in terms of the Fourier transforms 

^M^^.4^)-^-fe^^ /'°^7^,„ (7.19) 

1 - ^ll{u) 2 cosh ^ smh ^ 

(compare with eqn. (|7.7|)) and 

,-11 ^^^(0) -^-1 

1-$ll(0) t 

It coincides precisely with the corresponding limit S ^ 00 of eq. (|7.8| ) provided 

t - 1 (7.20) 

M = m 

We pause to discuss the logic. In the UV limit A = and we have a free boson after 
the shift dip dip + A. We can perturb this fixed point by A real or A purely imaginary. 
The first case is the usual sine-Gordon in the regime where there are only solitons in the 
scattering theory. The theory is massive, so the IR fixed point is the trivial one. The 
second case is like the sine-Gordon model with imaginary A, where there are only L and R 
"solitons" . As discussed in previous sections the theory is massless and the IR fixed point 
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is now nontrivial. Using known results about flow between minimal models (for example 
the large-t expansion discussed in sect. 2), we expect the IR fixed point to exhibit the shift 
t —>■ t — 1. In the IR limit, the LR scattering becomes negligible, hence the form of LL 
and RR scattering above. For LR scattering we require that in the opposite UV limit, 
the two scattering theories look the same, which is natural since they hold in two regimes 
connected at the UV fixed point. 

The expansion (7^) is obtained by analyzing our two sets of background field energy 



equations using a generalized Weiner-Hopf technique |]46|j30|] . This has been described in 



detail in so we do not present the full calculation here. Instead, we will explain how 
to extract the relevant information from the kernels. The technique relies on the usual 
Weiner-Hopf trick of dividing the Fourier transforms of the kernels into a product of two 
pieces, the first of which has no poles or zeroes in the lower half plane and the second none 
in the upper half plane. Defining l/K+{u!)K-{uj) = 1 — $(c<;), one finds that expressions 



of the form 



occur regularly in the analysis; the contour covers the upper half plane. The function 
h{uj) is different depending on where we are in the analysis, but it is analytic in the upper 
half plane. First we discuss the massive case. The poles in the contour are aX uj = i 
and at the zeros of 1 — which are at = 2m/(t + 1). The pole aX uj = i results in 
the bulk contribution £^(0). Ignoring the bulk term, (17^211) can be written as a series in 
exp(— 4i?/(t + 1)). In particular, the boundary condition results in an equation 



^e^ = const + V /i.5ne-^"^/(*+i) 



where and Qn are the residues of h{uj) / (cu— z)^ and g{uj), respectively. The hn themselves 
also obey an equation of this form, so for large A/M, we can write and each as a 
series in (A/M)"^/*^*"*"^). The energy is also given by a term like ( |7.21| ), so it too must be 
a series in (A/M )-'^/(*+^) as in (0): 

42 ^ /yi^N 4n/(t+l) 

£^^(AM)=^(«^)(0) Y.^n[^] . (7.22) 

^ n=0 ^ ^ 
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This gives the result (|7.4|) for ordinary sine-Gordon. For the imaginary couphng, we 
must first rewrite the equations ( [7.14]) in Weiner-Hopf form. The general result is that 

1 - $ll(0) 



The and obtained for the massless case are exactly the same as the ones ob- 
tained above for the massive one. The extra piece in the above expression for g{uj) is 
— sinh(7ra;/2)/ sinh(7rta;/2) here. It results in no other additional poles in the contour be- 
cause of the zeros in g. Its only effect is to change Qn to {—l)^gn (again ignoring the bulk 
piece). Since the above are not changed, we then find that the series for the massless 
flow is exactly the same as in ordinary sine- Gordon, except that the signs of every other 
term are different: 

4n/(l+t) 



with precisely the same coefficients /c„ as in expansion ( [7.22| ). 

We conclude that up to the known bulk vacuum energy contributions the massless BA 
system ( |7.14p -( [7.16| ) gives the correct analytic continuation of the massive one ([777|)-( [7r9D 
to purely imaginary ^, providing ( |7.20| ) holds. In particular, the low-temperature mass 
scale M is equal to the high-temperature scale m. 

This entire discussion can presumably be put on firmer ground by solving the Thirring 
model with imaginary mass using the traditional Bethe ansatz approach ME . 



8. Conclusions 

The examples presented here have mainly been the sine-Gordon and Thirring models, 
with and without a background field. These calculations can be extended to the truncated 
(RSOS) cases, the latter situation being completely unitary. See the references |lT^ , pl |. 



It seems that the physics of massless fiows is intimately related with the one of sym- 
metry breaking and that the conformally-invariant IR fixed points are generally some sort 
of Goldstone phase. Two of the simplest cases, the fiow from tricritical to critical Ising and 
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the flow from dilute to dense polymers have to do respectively with = 1 and N = 2 spon- 
taneous SUSY breaking (the latter being possible because of the non-unitarity) . Moreover, 
as already mentioned, the sine-Gordon model with imaginary coupling describes the flow 
from critical to low temperature 0{n) model (with n = 2 cos j). The Mermin- Wagner the- 
orem preventing spontaneous breaking of continuous symmetry does not apply to the cases 
n non-integer. As a consequence, 0{n) models in two dimensions —2 < n < 2 have a low- 
temperature phase which is massless and has properties reminiscent of a Goldstone phase 
||51j| (for instance they qualitatively agree with what can be deduced from the e-expansion 
in higher dimensions, extended formally to D = 2). 

Some important questions remain to be addressed. For instance, can one seriously 
describe conformal fleld theories using massless scattering (reconstruct Green functions 
using form factors)? As explained above, massless scattering is a sort of perturbation of 
the IR flxed point. Can one, using it (i.e. probably using the conserved quantities) give 
a meaning to the conformal perturbation theory of an IR flxed point by an irrelevant 
operator? 
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Figure Captions 



Figure 1: The structure of the ground state in the massive Thirring model. Left and right 
massless excitations are observed in the hmit ^ — > ±00. For instance, a pair of left and 
right solitons is obtained by making two holes for |^| ^ 1. 

Figure 2: Schematic structure of the ground state for imaginary uiq. 
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